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We study the time evolution of the entanglement negativity after a local quantum quench in (1+1 (-dimensional 
conformal field theories (CFTs), which we introduce by suddenly joining two initially decoupled CFTs at their 
endpoints. We calculate the negativity evolution for both adjacent intervals and disjoint intervals explicitly. For 
two adjacent intervals, the entanglement negativity grows logarithmically in time right after the quench. After 
developing a plateau-like feature, the entanglement negativity drops to the ground-state value. For the case 
of two spatially separated intervals, a light-cone behavior is observed in the negativity evolution; in addition, a 
long-range entanglement, which is independent of the distance between two intervals, can be created. Our results 
agree with the heuristic picture that quasiparticles, which carry entanglement, are emitted from the joining point 
and propagate freely through the system. Our analytical results are confirmed by numerical calculations based 
on a critical harmonic chain. 


I. INTRODUCTION 


respect to ,4-2’s degrees of freedom is defined as 


A. Introduction 


Recently, it has been recognized that quantum entangle¬ 
ment provides us a powerful tool to study quantum proper¬ 
ties of many-body systems in condensed matter physics m- 
3). When the system is prepared in a pure state |SR), a good 
quantity that describes the bipartite entanglement is the von 
Neumann entropy, which is defined as 

Sa = -Trp A lnp A (1) 

where pa = Tr bP is the reduced density matrix of subsys¬ 
tem A, with p = |q/) (SR |. An alternative measure of bipartite 
entanglement in pure states is the Renyi entropy 

= ^—^Trp n A . (2) 

1 — n 

These entanglement measures have proved to be of great use 
in characterizing quantum entanglement of many-body states. 

However, for a mixed state, neither the von Neumann en¬ 
tropy nor the Renyi entropy is a good measure of entangle¬ 
ment since quantum and classical correlations are not clearly 
separated in these measures. Now suppose we are interested in 
the entanglement between two subsystems A\ and A 2 , which 
are not necessarily complementary to each other and, are em¬ 
bedded in a larger system, the union Ai U A 2 cannot be de¬ 
scribed by a pure state after integrating out degrees of freedom 
in the complement of A\ U A 2 . In this case, we need to search 
for other quantities that may characterize quantum entangle¬ 
ment for a general mixed state. Among different proposals 
EH, a computable measurement of entanglement, the loga¬ 
rithmic negativity 0, turns out to be very useful and practi¬ 
cal. In particular, it is proved that the logarithmic negativity 
is a proper entanglement monotone in Ref. El- Following 
Ref. 0 , the negativity can be obtained by first taking a par¬ 
tial transposition. To be more precise, given a density matrix 
PAiUA 2 which describes a bipartite mixed state in a Hilbert 
space PLa^Ai = PLa x ® PLa 2 , the partial transposition with 
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where |e^) and \ef^) are arbitrary bases in PLa x and PLa 2 , 
respectively. Then the logarithmic negativity is defined as 


£ a u a 2 ■■= hr | |p 


t 2 

AiUA 2 


= InTijp 


t 2 I 

AiUA 2 I’ 


(4) 


where '/ 2 indicates the partial transposition with respect to A 2 , 
and the trace norm | l/TuuAn 11 is defined as the sum of all the 
absolute values of the eigenvalues of p l f LJ A ; . Recently, the 
logarithmic negativity has been extensively used to study vari¬ 
ous many-body systems, including one-dimensional harmonic 
chains EMQ), quantum spin chains EMI , free fermion sys¬ 
tems ED, and topologically ordered systems H9ll20l l. In par¬ 
ticular, the universal features of the entanglement negativity 
in one-dimensional critical systems have been understood by 
developing a conformal field theory (CFT) approach E1E2. 
Later on, the comparison of CFT results and numerical calcu¬ 
lations of one-dimensional critical systems were studied in a 
series of works j23ti25l . 

Although many works have been done on the entanglement 
negativity, there is less understanding on the non-equilibrium 
properties of the entanglement negativity. Most recently, time 
evolution of the logarithmic negativity after a global quench 
was studied with CFT approach l26l . In Ref. l27l . the nega¬ 
tivity evolution for two adjacent intervals after a local quench 
was numerically studied in a harmonic chain. However, a 
thorough study of the negativity evolution after a local quan¬ 
tum quench is still lacking, and it is appealing to unveil the 
universal features of the dynamical behavior of the entangle¬ 
ment negativity propagation. 

In this paper, our motivation is to study the time evolution 
of the entanglement negativity after a local quantum quench 
analytically. For simplicity, we consider a (l + l)-dimensional 
critical system, which is physically cut into two parts that are 
prepared in their own ground states. Then at time t = 0, we 
join the two parts together at their endpoints, and study the 
time-evolution of the entanglement negativity afterwards. As 
shown in Fig. [T] once the two CFTs are joined at the end¬ 
points, the interaction between them is introduced simultane- 
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FIG. 1: Setup for a local quantum quench. Two separate CFTs de¬ 
fined on two semi-infinite lines are joined together at their endpoints. 
Then quasiparticles, which may be viewed as entangled pairs, are 
generated at the jointing point, and propagate freely through the sys¬ 
tem. The entanglement negativity between two intervals which are 
far from each other may be built with the help of these propagating 
entangled pairs. 

ously, which generates quasiparticles (excitations) at the joint¬ 
ing point. These quasiparticles may be viewed as entangled 
pairs Il28]j32l which carry entanglement information. When 
the entangled pairs arrive at two intervals separately, the en¬ 
tanglement negativity can be built immediately. Because the 
(1+1) dimensional critical system is Lorentz invariant at the 
low energy limit, we can utilize the power of conformal field 
theory and understand the universal feature of this dynamical 
phenomenon. 

The rest of the paper is organized as follows. In part B of 
this section, we give a brief review of the path integral rep¬ 
resentation of the entanglement negativity, and then introduce 
the CFT setup for a local quantum quench in part C. In Section 
II, by using CFT approach, we compute the time evolution of 
the entanglement negativity for two adjacent intervals in part 
A, and two disjoint intervals in part B. We consider both cases 
where the two intervals are symmetrically and asymmetrically 
located. In section III, we describe the numerical method of 
calculating the entanglement negativity for a harmonic chain, 
based on which we study the local quench of the entangle¬ 
ment negativity. Then we compare the numerical results with 
the CFT results. In section IV, we conclude our work and list 
some interesting future problems to be studied. 

B. Entanglement negativity in quantum field theory 


interval r £ [0, /3]: 

P=7F f [dHx,T)}Y[5{(l>{x,0) - 

x ( 6 ) 

x 0<5(#r,/?) - (/>"{x"))e Se , 

X 

where the rows and columns of the density matrix are labeled 
by the fields {</>(£, r)} at r = 0 and /3 respectively, with /3 
being the inverse temperature, Se is the Euclidean action and 
Z = Tre >n is the partition function. Now we consider sub¬ 
systems A 1 and A 2 located in intervals [u^Ui] and [u 2 ,v 2 ], 
respectively. Then the reduced density matrix p Al uA 2 may be 
obtained by sewing together all the points along edges r = 0 
and t = j3 except the points in A\ U A 2 . That is, we leave two 
open cuts at [ui, iq] and [u 2 , v 2 ] along r = 0. 

Next, before we compute Tr (pJiiM,)™'’ ^ * s beneficial to 
see how to calculate Tr ( PA 1 uA 2 ) n first- In order to calculate 
Tr ( pA 1 uA 2 ) n , we consider n copies of the cut plane, and then 
sew together the cut [it*, Vi] 3 T=0 _ with the cut [u^, t’i]^ 0+ for 
* = 1,2 and all the copies j = 1, • • • , n. Note that for j = n, 
we sew together the cut [it*, q_ with the cut [m, Uj]^Z 0 + - 
In this way, we define a n-sheeted Riemann surface TZ n . The 
trace of {pA 1 uA 2 ) n is then given by 

Tr ( PAiUA2 ) n = ^, (7) 

where Zjz n is the partition function for the orbifold CFT on 
lZ n . Rather than dealing with the fields on a nontrivial mani¬ 
fold, it is found more convenient to work on a single complex 
plane. It turns out Eq. 0 can be expressed in terms of lo¬ 
cal twisted fields defined at (ui, 0) and (vi, 0) on the complex 
plane as follows 

Tr ( PA 1 UA 2 ) n = ('T n (ui)Tn(vi)Tn(u 2 )T n (v2 )} ■ (8) 

Intuitively, the effect of twist fields T n and T n is shown in 
Fig-0 Winding anticlockwise (clockwise) around the twist 
field T n (Tn), once the branch cut is crossed, one will go from 
layer j to layer j + 1. 

With the introduction of twist fields, the expression of 
Tr (fft uAo ) n is very straightforward. As discussed in Refs. 
OT1 l22l . the effect of partial transposition with respect to A 2 
is equivalent to changing the two twist operators T n (u 2 ) and 
T n {v 2 ). Then one has 


A detailed description of path integral representation of the 
entanglement negativity can be found in Ref. Il22l . For the 
completeness of this paper, we give a brief review here. 

First, as discussed in Ref. Il22l . by using a replica trick, one 
can relate the entanglement negativity with the integer powers 

of Pa\ua 2 as 


£ Ai ,a 2 = f lim i lnTr(pJ iUA2 ) , (5) 

where n e is an even integer, and the density matrix p may be 
expressed as a (Euclidean) path integral in the imaginary time 



(Tn(ui)Tn{Vl)Tn('U2)T n (v 2 )) . 


(9) 


If the two intervals [iti,Ui] and [ 142 ,^ 2 ] ar e adjacent to each 
other, we simply set u 2 —> tq , and then Eq. ([9]) can be written 
as 


Tr 


(PMUAiY = (Tn(ui)fZ(u 2 )Tn(v 2 )) . (10) 


Therefore, from Eqs. ([5]), |9| and (10 1 , it is found that the com¬ 
putation of the entanglement negativity reduces to the compu¬ 
tation of expectation values of twist fields in a complex plane. 
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FIG. 2: Path integral representation of (a) Tr(pA)" and (b) Tr(pJ) 71 for two disjoint intervals. 


C. CFT approach to a local quench 

Before we study the CFT approach to a local quantum 
quench, it is beneficial to comment on the difference between 
local quenches and global quenches. Local quenches are more 
complicated than global quenches because they are inhomo¬ 
geneous. For global quenches, we change the parameters of a 
translational invariant Hamiltonian globally, and therefore the 
system before and after global quenches are always transla¬ 
tional invariant. In this case, as discussed in Refs. 01113211 , the 
initial state can flow to a conformal invariant boundary state 
under renormalization group (RG). For local quenches, how¬ 
ever, before we join the two decoupled CFTs together, the to¬ 
tal system is apparently not translational invariant. Therefore, 
the initial state cannot flow to a conformal invariant boundary 
state. In addition, for global quenches, quasiparticle excita¬ 
tions are emitted from everywhere in the bulk of CFT; for lo¬ 
cal quenches, quasiparticle excitations are emitted only from 
the point where two CFT are joined together. 

The time dependent density matrix can be written as p(t) = 
\^{x,t)){cj>{x,t)\, where | (j>(x,t)} = e~ lHt \(po(x)). In path 
integral representation, one has 

where the factor e~ eH is introduced to damp out high-energy 
modes and make the path integral absolutely convergent. If 
the CFT arises as the low energy limit of a lattice model, then 
e may be viewed as the lattice spacing. In the study of global 
quenches 0X14331 . \<j)o(x)) may be considered as a conformal 
invariant boundary state under the RG. For local quenches, 
as we discussed above, the initial state cannot flow to a con¬ 
formal invariant boundary state under the RG. In this case, 
one may introduce boundary condition changing operators, as 
utilized in Refs. [34j :351. In this work, however, we will fol¬ 
low the method proposed by Calabrese and Cardy 0X1 . As 
shown in Fig[3] the density matrix can be expressed in terms 


of the path integral on a modified word-sheet, where the phys¬ 
ical cut corresponds to two ‘walls’ with one extending from 
r = —oo to —e and the other extending from r = +e to +oc 
in a complex z-plane. No energy nor momentum can flow 
through the two ‘walls’, and therefore conformal boundary 
conditions are imposed on the wall (As will be shown later, 
the concrete boundary condition does not affect the universal 
result we consider.). For convenience of calculation, we map 
the z-plane to a right half plane (RHP) with Re w > 0 by 
using the conformal mapping 


w = 



( 11 ) 


Then the local quench problem is reduced to the calculation 
of correlation functions of twist fields in the RHP 0X1 . which 
we will study in detail in the next section. 


2s 
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FIG. 3: Illustration of the conformal mapping in Eq.(|TTJ, based on 
which the z-plane is mapped to a right half plane (RHP) with Rew > 
0. For later use, we also label Wi = —u>i, which is the image of wi. 
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FIG. 4: Configurations of two intervals considered in this work: two 
adjacent intervals (up) and two disjoint intervals (bottom). 


II. ENTANGLEMENT NEGATIVITY AFTER A LOCAL 
QUENCH: CONFORMAL FIELD THEORY APPROACH 

In this section, we calculate the time evolution of the en¬ 
tanglement negativity after a local quench in conformal field 
theories. We will consider adjacent intervals in part A and 
disjoint intervals in part B, respectively. 


A. Two adjacent intervals 

1. Semi-infinite intervals 

As a warm up, we consider the simplest case, i.e., the total 
system is bipartitioned into two semi-infinite parts A \ and A ■>. 
In this case, pa 1 ua 2 is pure, and the logarithmic negativity is 
the same as the Renyi entropy with n = 1/2 1 21 1 [22j. This 
case was also studied in Ref. ED- 

For two adjacent semi-infinite intervals, we only need to 
consider a single twist field 7^(zi) in z> -plane, which is in¬ 
serted at 


Z\ = l + IT. 


( 12 ) 


By choosing the insertion position at the origin l = 0, i.e., 
A i £ (—oo,0] and A 2 £ (0, +oo), one simply has z\ = ir. 
The expectation value of 7^(zi) can be expressed as | 


<T„ 2 ( 2 i)) = c n 


dw 


dz 


A (2) 


2Re(wi) 


(13) 


where c n is a nonuniversal constant which depends on the par¬ 
ticular boundary CFT, a is an UV cutoff (e.g., the lattice spac- 
ing in a harmonic chain), and A„ ' is the sca ling dimension of 


Tn{z). By using the conformal map in Eq. (Ill, one has 


Wi 


= i- + ~ \/e 2 - 


(14) 


and 


dw 


v/e 2 — • 


With an analytical continuation r —► it, one can obtain 

A< 2 ) 






2 (e 2 + f 2 ) y 


(16) 


As discussed in Refs. ED EH, the scaling dimension A], 2 ' 1 
depends on the parity of n as 


Ai 2 ) = 


An 


odd n. 


2 A n /2 even n, ’ 


where 


A c ( 1 

A n = -r n - - 

12 V n 


(17) 


(18) 


Then by using the expressions in Eqs. ([5]) and ([9]), one can get 


£ = — In 


ae 


4 2(e 2 + 1 2 ) 


+ c lt 


(19) 


where c] = lnci. As in Refs. 01331, the short time behavior 
of £(t) allows us to fix the regulator e in terms of the non¬ 
universal constant ci by requiring 


f(f = 0) = -|ln^ + c;=0, 


based on which one gets 


e = -e _ 4 a i/ c 
2 


( 20 ) 


( 21 ) 


Now it is possible to eliminate a and c\ in Eq. (19 1 in terms of 
e, and then one can get 


£ 4 ln e 2 + t 2 ' 


In the limit t e, one ends with 

£ C 1 1 

£ =2 ln e 


( 22 ) 


(23) 


which was observed in the numerical calculations based on a 
critical harmonic chain ED- 


2. Symmetric finite intervals 

In this part, we consider the case of symmetric finite inter¬ 
vals with l\ = l 2 = l, i.e., Ai £ [—1, 0] and A 2 £ (0, l ], as 
shown in Fig. [4] In this case, pa, ija 2 represents a mixed state, 
and there is no correspondence between the logarithmic neg¬ 
ativity and the Renyi entropies. By using Eq. ( flQ| and doing 
a conformal mapping onto the RHR one has 


Tr 


(P5U)" = (Tn(Zl)TZ(Z2)T n (z 3 )) 


n 

i=l 


dw 

dz 


Hi) 


(Tn(w 1 )TZ(w 2 )T n (w 3 )\ 


dz 


( 15 ) 


(24) 
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FIG. 5: Entanglement negativity £ for two symmetric adjacent intervals as a function of time. Here we choose the central charge c = 1, 
e = 0.1, l = 25, 50, 75 and 100, respectively. Shown in (a) is the CFT result, and (b) is the numerical calculation based on a critical harmonic 
chain. 


where the scaling dimensions Ah) = A( 3 ) = A n and A( 2 ) = 

A), . The three-point correlation function on the RHP can be 
expressed as 


(Tn{wi)Tn(w 2 )Tn(w 3 )) 


RHP 

A (2 >-2A„ 

Vl,3 


1/2 


n;=i ik -Wi)/a 


A 


«) 


A 


( 2 ) 


A 


( 2 ) 


v 7l,2 ^2,3 


•F{{Vj,k })> 

(25) 


where rjij are cross ratios which can be constructed from the 
endpoints Wi (and their images w L ) of the intervals in the RHP 
as follows 


which may be further expressed as 


e = 



e 2 +1 2 


C W13W12W23 

4 W^W 12 W23 ’ 


(28) 


where we have defined Wij = \wt — Wj |, w^ = \wi~Wj \ and 
w ij = \wi — Wj\, respectively. 

With the expressions of that are calculated in the Ap¬ 
pendix, one can obtain the entanglement negativity £ as a 
function of time as follows 


£ = 




l -f -1 
~t ’ 




t < l 
t > l. 


(29) 


_ _ {Wi-Wj){Wi-Wj) 

'h,j — ( ~ \ 

[Wi - Wj ) \Wi - Wj) 

with Wi = —Wi being the image of (see Fig|3]>. The 
nonuniversal function F({r/j.fc}) depends on the full operator 
content of the CFT. J r ({ r lj ! k}) is usually difficult to calculate 
and only known for several specific CFTs and BCFTs. But 
it is found that in the limits 77 ^ —► 0, 1, or 00 , the function 
J-({ilj t k}) is just a constant, which follows from the long- and 
short-distance expansions of the correlation functions of twist 
operators |[3l [32] 13614381 . For symmetric intervals in this part, 
we calculate the cross ratios r/ij explicitly in the appendix. It 
is found that one always has rj l:) = 1 or 0 for the cases fCl, 
t = l + 0~ and t > l. In other words, our results are universal 
for the above three cases. 

By using Eqs. (|5]i and ( flQ| ), and neglecting various nonuni¬ 
versal terms, we have 


E 


C , £ 2 C T? 1,3 

- In —---In--, 

4 e 2 + 1 2 8 t?i 2^72,3 


In the limit l, t /§> e, £ can be simplified as 


£ = 


c t 2 + e 2 (l -t) 
4 ~(K)’ 

c l 

— in — + const, 


t < l 
t > l. 


(30) 


Shown in Fig. [5ja) is the plot of E with different l. At the very 
beginning of the local quench t <C l, based on Eq. (30 1 , one 
has 


C C t 

E= 2 ln ? 


(31) 


which agrees with the result of semi-inhnite intervals as 
shown in Eq. ( |23] l. This is reasonable because in the limit 
t < /, the quasiparticles essentially propagate without notic¬ 
ing the finite size effect. For t < l, the entanglement saturates 
for a certain time. Then for t > l, we get the ground-state 
value of £, i.e., 


£a = jK 


(27) 


(32) 




























6 


which is also observed in the numerical calculations in Ref. 
03. Note that in the numerical calculations, £ tends to the 
ground-state value gradually. In our CFT results, £ drops 
to the ground-state value immediately after t = l. This is 
because all quasiparticles propagate at the same velocity in 
CFTs. In lattice models, however, the dispersion relation is 
not linear for all momentum vectors, and therefore not all 
quasiparticles propagate at the same velocity, as discussed in 
detail in section IV. 

In addition, the scaling behavior of £ for a harmonic chain 
was numerically studied in Ref. Il27l . For t < Z, they proposed 
the ansatz 


On the other hand, from the scaling behavior in Ref. Il27l . one 
can find 

Numerical (i = ^ ) = \ In ^ + 0.131n3 - 0.28 lu2. (35) 

In the large l limit, i. e ., l e, one always has 

£cFT (t = —) ~ £numerical (f = ~ In (36) 


£ ~ In 


t a {l-t)P 
l~P (l + t 


(33) 


By fitting the numerical results, they found a = 1/2, (3 ~ 
0.15, 7 ~ 0.13 and p = -(/? + 7 ) ~ -0.28. For our CFT 
results in Eq. (|30j», by setting c = 1 and taking the limit f > e, 
one has a = 1/2, (3 = 0.25, 7 = —0.25 and p = — (/3 + 
7) = 0. On the other hand, in the limit t -C L both Eq. 
(30 1 and Eq. (33 1 collapse to £ = / In t. We attribute the 
above disagreement/agreement to the following fact. For t < 
l, because we neglect the non-universal functions J r ({r / 7 fe }) 
which may not be constants, our results are not accurate and 
therefore may not obtain the correct scaling behavior. For t <C 
l, however, our CFT results are universal and independent of 
the specific CFT. To reproduce the numerical results for t < l 
in Ref. E3- we have to consider the nonuniversal functions 
^{{ r lj,k})> which is a difficult task, and out of the scope of 
our work. 

Nevertheless, by comparing the values of the plateau be¬ 
tween CFT results and numerical results in Fig. [5] it is found 
they are very close to each other. To be concrete, let’s take 


£(t = for example. Based on Eq. (301, one can get 


cm 


— In — 

+ const. 


2 2e 


- - In 

e 2 

-In 

e 2 + 1 2 

4 

8 

- - In 

e 2 

-In 

e 2 + 1 2 

4 

4 

- - In 

e 2 

- In 

SO 

+ 

(M 

4 

4 


(34) 


Before we end this part, we mention that it is interesting to 
check how £{t) behaves in other lattice models. Considering 
£{t < l) depends on the non-universal functions 
which varies for different CFTs, we expect that for other crit¬ 
ical lattice models such as the critical Ising model one may 
observe different scaling behaviors in £{t < l). 


3. Asymmetric finite intervals 

In this part, we consider the case of asymmetric finite inter¬ 
vals with A\ £ [— 1\, 0] and A 2 £ (0, 12 ], as shown in Fig. [4] 
Without loss of generality, we suppose l\ < 1-2 . The calcula¬ 
tions are similar to the symmetric case, and we need to eval¬ 
uate the three point correlation functions in Eq. First, as 
shown in the appendix, we calculate the cross ratio rjij explic¬ 
itly. It is found that one always has rjij = 1 or 0 for the cases 
t <C Zi, t = I 2 + 0~ and t > ( 2 - That is to say, our results 
are universal in these regions. Second, by neglecting various 


nonuniversal terms, we arrive at the same result as in Eq. (28 1 . 


The difference is that for asymmetric intervals, we have dif¬ 
ferent expressions of w t;] , as explicitly given in the appendix. 
By plugging w,j :j into Eq. (28 1 , one obtains 




eh l 


t < l[ 
h < t < I 2 
t > ( 2 * 


172 


In the limit l,t e, £ can be simplified as 


-In 

(■h - t){l 2 - t)t A 

t < h 

8 

(h + t)(l 2 + t)e A 


hi 

( h — h)l\t 2 

A 

c-+- 

A 

8 

(l\ + h ){t + h){t — h) 


7 ln 

hh 

— -— + const, 

t > / 2 - 

4 

e(h + l 2 ) 



Note that in the limit h = I 2 , we can reproduce the symmetric 
adjacent intervals result in Eq. pO} . The plot of £ for different 
{h,h) is shown in Fig.^a). First, for t <C min[Zi,? 2 ]> one 
can find the time evolution of £ (t) is similar to the symmetric 
case. Actually, based on Eq. one can check that in the 
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FIG. 6: Entanglement negativity £ for two asymmetric adjacent intervals as a function of time. Here we choose central charge c = 1, e = 0.1, 
(h, h) = (25, 75), (50, 75), and (75, 75), respectively. Shown in (a) is the CFT result, and (b) is the numerical calculation based on a critical 
harmonic chain. 


limit t -C min^!, l 2 ] 

£ = |ln~, (38) 

which shows the In t behavior again, as expected. For t > 
max[Zi, l 2 \, we obtain the ground state value of the entangle¬ 
ment negativity EDEU, i.e.. 


_ c hi 

~ 4 n e(h + l 2 )' 


(39) 


Interestingly, it is found that the sudden drop of £ happens at 


t = min[Zi, l 2 ], 


(40) 


which is again straightforward to understand based on the 
heuristic physical picture that the quasiparticles may be 
viewed as entanglement pairs of two quanta. For t < 
min[Zi, l 2 ], two entangled quasiparticles are in A\ and A 2 sep¬ 
arately, and create the entanglement between A\ and A 2 . At 
t = min[Zi, l 2 ] = h (here we suppose Zi < l 2 ), although one 
quasiparticle is still in A 2 , the other quasiparticle propagates 
out of A 1 , and therefore the entanglement between A\ and A 2 
decreases suddenly at t = inin[/-|, l 2 \. 

Before we end this part, we emphasize that our results are 
universal for the regimes t. <C h, t = l 2 + 0“ and t > l 2 . For 
the case of t. < l 2 , however, similar with the symmetric case, 
because we neglected the nonuniversal functions ^({Vj^}) 
which may not be constants, our results are not accurate, and 
one has to calculate ^({ry^ fc }) for different CFTs. 


B. Two disjoint intervals 

1. Symmetric finite intervals 

In this part, we consider the symmetric disjoint intervals, 
i.e., A\ £ [—d — l , —d] and A 2 £ [d,d + l ], as shown Fig. 


|4]b). In this case, we need to consider the correlation func¬ 
tion of four twist fields as shown in Eq. & By applying the 
conformal map in Eq. (fTT), one has 


{Tn^Tniz^Tniz^TniZi)) 


4 dw A " 


n 


dz 


{Tniw^Tniw^Tniw^TniWi))^ 


where the four-point correlation function on the RHP has the 
form 

(T n (Wi )T n ( W 2 )T n (w 3 )T n (u> 4 ) ) RHp 

-2 1 / \ A ( l 2 ) / 2 -A„ 

= _<_ 1 / 771 , 4772,3 \ 

nil \( w i - W,;)/a| A ™ Vl, 2 V 3,4 \Vl,3V2,4,J 
X ({Vj,k}) ■ 

(41) 


For the nonuniversal functions ^j), as explicitly cal¬ 

culated in Ref. l36l . they are simply a constant in the limit 
l/d <C 1. In other words, when the two intervals are far apart, 
we do not need the knowledge of F ({%*;})• By using the 
definition in Eq. 0- and dropping various multiplicative con¬ 
stants, we have 


£ = -®l 

8 \ ^1,3*72,4/ 

which is alternatively written as 


(42) 


C W 1 ±W 1 2W23W 2 l W ii W i3 W 23 W 2i 
8 W14W13 ^23^24^14^13^23^24 ' 


(43) 


By noting that Wij 
simplified as 


w ij and w ij = 


w---, Eq. (431 can be 


C Wi4'U> 1 3W 2 3'U>24 

4 W 1 4’U)13'IU 2 3'UI24 ’ 


( 44 ) 
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Time Time 


FIG. 7: Entanglement negativity £ for two symmetric disjoint intervals as a function of time. Here we choose the central charge c = 1, t = 1, 
(d,l) = (140,10), (160,10) and (180,10), respectively. Shown in (a) is the CFT result, and (b) is the numerical calculation based on a 
critical harmonic chain. 


With the explicit forms of Wij given in the Appendix, we can 
obtain the entanglement negativity £ as a function of time as 
follows 


0 


£ = 


c (2 d + l)(d + l — t)(t 2 — d 2 ) 
4 edl(d T l + t ) 

c (2 d + l) 2 

4 4 d(d - 


0 


t < d 


d <t < d + l 


t + d + l 


(45) 

Note that in the study of the negativity evolution after a global 
quench, it was found that £(t) shows the same behavior as 
the Renyi mutual information apart from the prefactor (26). 
For the local quench studied here, by comparing our result in 
Eq. ( [45) with the result of mutual information in Ref. (36), it 
is found that the expressions are also the same except for the 
prefactor. In other words, our results parallel with the story in 
the negativity evolution after a global quench. The relation be¬ 
tween the entanglement negativity and the mutual information 
after a local quantum quench will be systematically discussed 
in section IV. 

As shown in Fig. [7) a), we plot the evolution of the entan¬ 
glement negativity with different {d, l) according to Eq. (45 1 . 
A Tight-cone’ effect can be observed: For t < d, there is 
no entanglement negativity between A\ and A 2 - At t = d, 
the entanglement negativity begins to develop, and reaches the 
maximum approximately at i = d + l/ 2. At t = d + l, the en¬ 
tanglement negativity decreases suddenly, which corresponds 
to the entangled pairs propagating out of intervals A t and A 2 


simultaneously. Note that at t = d + 1/2, taking the limit 
d^$> l, one has 


^t=d +1 



(46) 


which is independent of the distance d, as also can be ob¬ 
served in Fig. [7] That is to say, with the help of entangled 
pairs, we can create a long-range entanglement between two 
intervals which are far from each other. Note that this long- 
range entanglement was also observed in the time evolution 
of mutual information I(t) in Ref. (36), where it is found that 

T ^ C 1 l 

1 t=d+ 4 “ 3 111 2e ' 


2. Asymmetric finite intervals 

In this part, we consider the asymmetric disjoint intervals. 
We have multiple choices as follows: (i) d\ 7 ^ d 2 , h = h, (ii) 
di = d 2 , h 7 ^ 1 2 and (iii) d\ 7 ^ d 2 , l\ 7 ^ l 2 . For simplicity, 
we consider the case (i), i.e., A\ £ [— d\ + l , —di] and A 2 £ 
[d 2 , d 2 + Z], Without loss of generality, we choose d\ < d 2 < 
d\ + l. 

The calculation of the negativity evolution is similar with 
the symmetric case, and we obtain the same result in Eq. (44) . 
The difference is that we should express w^j in terms of g^, 
d 2 and l, as explicitly shown in the appendix. By plugging 
the expressions of w t j into Eq. (44) , one arrives at the time 
evolution of the entanglement negativity 































9 




FIG. 8: Entanglement negativity £ for two asymmetric disjoint intervals as a function of time. Here we choose the central charge c = 1, e = 1, 
l = 15, (di,d, 2 ) = (150,150), (150,155), (150,160), and (150,165), respectively. Shown in (a) is the CFT result, and (b) is the numerical 
calculation based on a critical harmonic chain. 


0 


£ = 


c (di + — t)(c ?2 + l + f )(<^2 — d\ + V) 

8 (di + c?2 + 0(^2 + l — t) (d‘2 + i)(d2 — d\) 

c e(di + (^ 2 ) / (d\ + l — X 2 XX 2 T Z — di)(di -(-1 + t)(d 2 + l + t) 

4 n 2(di + d 2 + l) y (di + Z - f)(d 2 + l ~ t){t 2 - d\)(t 2 - d^) 

c [f 2 - (di + Z) 2 ](di + d 2 ) 2 (d 2 + l — di)(di + d 2 + 21) 

8 (d 2 - di)(di + d 2 + l) 3 (t 2 — d 2 ) 

c (di + d 2 + 2Z)(di + d 2 ) 

4 (di + d 2 + l) 2 


t < di 
di <t < d 2 

d2 < t < di + l 

di4-Z<f<Cd2 _ t _ Z 

t > d2 4“ l 


(47) 


One can check that when di = d 2 = d, the result in Eq. 
is reproduced. 

According to Eq. (47 1 , we plot £(Z) with different (di, d 2 ) 
in Fig. [8})a). Compared to the symmetric case, the ‘light-cone’ 
effect is still observed. The difference is that the time when 
£ (t) increases quickly now happens at 

t = maxfdr, d 2 ], (48) 

and the time when £ (t) decreases quickly happens at 

t = minfdr + l, d 2 + l], (49) 

which is also in agreement with the quasiparticle picture. 


entanglement negativity for a harmonic chain has been nu¬ 
merically studied in several works lf9l l22l l26l l27l 39 1. Here 
we follow the method developed in these works, and apply it 
to the local quench problem. We will first introduce the lat¬ 
tice model and the covariance matrix in part A. In part B, we 
introduce the evolution matrix and show how to calculate the 
entanglement negativity. In part C, we apply the method to 
the cases studied with CFT approach, and compare the results 
accordingly. 


A. Harmonic chain and the covariance matrix 


HI. NUMERICAL EVALUATION OF THE NEGATIVITY 
FOR A HARMONIC CHAIN AFTER A LOCAL QUENCH 

In this section, to confirm our CFT results, we study the 
time evolution of the logarithmic negativity after a local quan¬ 
tum quench in a lattice model, a critical harmonic chain. The 


The Hamiltonian of the harmonic chain is 

H = Dll + + y (9n+1 - 9n) 2 ], (50) 

n— 1 

where N is the number of sites of the chain, M is the 
mass scale, wq is the characteristic frequency, and I\ is the 
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nearest-neighbor coupling. p n and q n denote the momentum 
and position operators with canonical commutation relations 
\Pm Pm] = \0_rn Qm] = 0 and [ q n , Pm] = '^n.m ■ 

For periodic boundary condition (PBC), the Fourier trans¬ 
form of the canonical variables are 


Qn 


=£ 


i 


'- qk V~L 


0 2ivikn / L 


(51) 


»=l>7r 


3 —27t ikn/L 


where n = 1, • • • , L. For p n , the Fourier transform is identi¬ 
cal to q n . The Hamiltonian is diagonalized in the momentum 
space 


( 1 ^2 , MuJ k ■$. 


^ = £ 


k= 1 


\2M 


Pk 


Qk 


(52) 


where 


4 K . ,nk , 


= \ful + — sin(—) 2 , A; = 1, • ■ ■ ,L, (PBC). 


M 


(53) 


For the Dirichlet boundary condition (DBC), the Fourier 
transform is not valid due to the breaking of translational sym¬ 
metry. However, the Fourier sine transform can be defined as 


L—l 

9™ = £ % 
k =1 
L—l 


2 ( 7 vkn\ 

qk \lL sm {~ir) 


Qk 


(54) 


E l 2 ( 7 T kn\ 

Qn.\l — sin I — 1 


\ l r 


where n = 1, ■ ■ • ,L. For p n , the Fourier sine transformation 
is defined similarly. The Hamiltonian in the momentum space 
is identical to Eq. ( |52| . But the frequency u> k has a different 
form. 


4 K 


, 7T k . 


u k = \/uj$+ — k = l,---,L-l, (DBC). 

(55) 

The covariance matrix is constructed from the two-point 
correlators 


_ t> p f (0\q n Qm\0) (0|(?nP m |0) \ 

7n, m - ne ^ (0 | pngm | 0) ( 0 |p n p m |0) ) ' 


(56) 


For PBC, the correlators are 


1 ^ ^ 

< 0 M™| 0 } = -]T — cos 


k =1 
L 


27r/c(n — m) 


( 0 |p n Pm| 0 ) = — ^ Mw fc cos 


fc=l 


L 

27rfc(n — to) 


(0|<ZnPm|0) — Zf) n7n /2. 


(57) 


For DBC, the correlators are 


L ^ Mm, Sm \ L ) ™ ( L ) 
(0\p n Pm\0) = ^ £ Mutfc sin sin , 

(0|g„p m |0) = 2. (58) 

B. Evolution matrix and the logarithmic negativity 

From the Heisenberg equation of motion, q k [t) = jjPk(t) 
and p k [t) = —MuJ k qk(t), we have 

Qk{t) = (cos u} k tq k ( 0 ) +w i : 1 sinu;fctp fe ( 0 )), 

Pk{t) = VM(-u) k sinu} k tq k (0) + cosoj k tp k (0)). 

The time-dependent canonical variables in the real space are 
1 


' V 

Qn(t) 2^ <t > t( n )<t>k{ m ) 


x (cosWfctg m ( 0 ) + uj k 1 sinu;fetp m ( 0 )) 
Pn(i) =v / M£^(n)0 fc (m) 


k,m 


x (-w fc sinwfet 9 m ( 0 ) + cosw fc fp m ( 0 )), 


where 


^n(fc) = -Le- 27 ^"^, 


(PBC), 


((>„(&) = —j= sin( 7 rfcn/L), (DBC). (59) 

\JL 

Therefore, the time evolution of the covariance matrix is 

7 (f) = 5'(t)7(0)5'(t) T , (60) 

where the evolution matrix is 

Sn„(t)=Y. 4>l(n)<t>k{rn) 

k 

x ( VM C0SUJkt 

l — yfMuJk sinwfct '/M cosu> k t I 


( 61 ) 


The entanglement properties are encoded in the reduced 
density matrix, which can be extracted from the the covariance 
matrix 74 associated with the the subsystem A. The logarith¬ 
mic negativity is defined by the partial transposition of the 
reduced density matrix pa with the subsystem A = Ai U A -2 
as £ = InTr | p A 2 |. We first consider the partial transposition 
of 74, which can be constructed by inverting the signs of the 
momenta corresponding to A 2 ]9l . 


1 T a = 


h A 0 !a 

0 i A ®U 2 


• 1A ■ 


h A 0i A 

0 Ia 


(62) 
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where 


1 0 
0 -1 


n,m € A 2 , (63) 


and l a is the number of sites in region A. The symplectic 
spectrum of 7 J 2 can be obtained by exact diagonalization after 
multiplying with a symplectic matrix E 


E = 


0 

-h A 



(64) 


with uj k = y Uq + 7 ^ sin(|^) 2 , and ”/i/ 2 (t) denotes the co- 
variance matrix of the disconnected chain. The evolution ma¬ 
trix in the situation is 


»(*) 



x 


777 cos Q k t V 7 vf^fc sinf 2 fc f\ 

sinflfcf \/Mcosflfcf J ’ 

(69) 


After computing the spectrum of zE • 7 a with the set of 
eigenvalues {±A 1; ±A 2 , • • • ± A i A ; A M > 0}, the trace norm 
of the partial transposition of the reduced density matrix is 

EM 


Ia 


Tripoli 


I'V' + 2 I _ 1^ _ 2 I 



Therefore, the logarithmic negativity is 


Ia 

£ = - In [min(l, 2A /J )]. ( 66 ) 

w=i 


where £l k = yju> 1+ sin(^) 2 . 

In the following, the negativity is computed by setting 
M = K = 1 and ojq = 0 so that the maximal group ve¬ 
locity of normal-mode excitations is set to unity. (Hereafter 
we will use the ‘light-speed’ v c (= 1 ) to represent the maxi¬ 
mal group velocity.) The total length of the harmonic chain 
is N = 500. The partial transposition is performed with re¬ 
spect to the interval A 2 . Notice that for u ;0 = 0, the system is 
critical with the central charge c = 1. As shown in Figs.[5]'-[8] 
we compute the negativity evolution for both adjacent and dis¬ 
joint intervals, including symmetric and asymmetric cases. By 
comparing the analytical results obtained from CFT approach 
and the numerical results based on the harmonic chain, it is 
the found that the main features agree very well. 


Let us summarize the details of computing the negativity: 


Constructing the covariance matrix based on Eq. (56 1 . 


D. Comparison between CFT results and numerical results 


Constructing the evolution matrix S rl rn (t) [Eq. (61 1 ] 
and obtain the time-evolution covariance matrix 7 (f) 
[Eq. @], 


• Constructing the reduced covariance matrix 7.4 (t) and 
computing the spectrum of zE • 7 . 4 (f). Then, the nega¬ 
tivity is obtained by Eq. (| 66 |. 


C. Evolution of entanglement negativity after a local quench 


We consider two disconnected harmonic chains with equal 
number of sites N under the Dirichlet boundary condition, 
with each harmonic chain prepared in its ground state. At 
t = 0 , two disconnected chains are joined together as one 
harmonic chain with number of sites 2 N under the Dirichlet 
boundary condition (similar to Ref. Il27l ). 


From Eqs. (56 1 and (58 1 , the covariance matrix at t = 0 is 


7(0) - ( 7l(0) 0 

7(Uj - l 0 72 ( 0 ) 


(67) 


where 


[7l(0)]n,m = [72 (0)]n,ri 


1 


N -1 


/ Trkn\ 

N - sm sin {-N 


=-Y 


( Trkm\ 


0 


1 

Mui k 

0 Mu} k 


( 68 ) 


Although the main features of our CFT results and numeri¬ 
cal results agree with each other, we notice that there are some 
disagreements in detail between the two methods as follows. 

(i) For the results of two adjacent intervals in Fig. [5] at time 
f = l/v c = l, the entanglement negativity obtained from the 
CFT approach drops to the ground-state value suddenly. For 
the lattice model, however, the entanglement negativity ap¬ 
proaches the ground-state value gradually. This is due to the 
fact that in CFTs all the quasiparticles propagate at the same 
speed v c , which is not the case in a lattice model. As shown in 
Eq. ( |55j ), the dispersion relation is nonlinear, which indicates 
that not all the quasiparticles have the same group velocity. In 
particular, for quasi-particles with higher energy, their group 
velocities v are smaller than the light speed v c , i.e., v < v c . 
As shown in Fig. [9] qualitatively, we can divide the entan¬ 
gled pairs into four groups according to their group velocities 
(■ vl,vr ) : 

fast-fast pair 
fast-slow pair 

(70) 

slow-fast pair 
slow-slow pair. 

In the CFT study of two adjacent intervals, there are only fast- 
fast pairs which lead to the abrupt drop of £ at t = l. In the 
lattice model, however, the slow-slow pairs still make contri¬ 
butions to £ even for t > l, and this is why £ drops to the 
ground-state value gradually in our numerical results. 


(' v L ,v R ) ^ 
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FIG. 9: Effects of propagating entangled pairs on the entanglement 
negativity £ of two adjacent intervals and two disjoint intervals in a 
lattice model. For two adjacent intervals, all the four kinds of en¬ 
tangled pairs may contribute to the entanglement negativity. For two 
disjoint intervals, however, only the fast-fast entangled pair may con¬ 
tribute to the entanglement negativity near t ~ d, as discussed in the 
main text. 


FIG. 10: Distribution of entangled pairs as a function of ( vl,vr) 
for a harmonic chain. The peak at (—vl/v c ,vr/v c ) — (1,1) is 
contributed by the fast-fast pairs. The region ( —vl/v c , vr/v c < 1 ) 
with finite £ is contributed by the other three kinds of entangled pairs, 
i.e., fast-slow pairs, slow-fast pairs and slow-slow pairs, which result 
from the high energy contribution. The parameters we choose are 
^ = l 2 = l — 5 and f = 100. 


(ii) The concrete values of the entanglement negativity ob¬ 
tained from the CFT approach and the numerical method do 
not agree with each other in a perfect way. In particular, 
the CFT results are much larger than the numerical results 
for the disjoint cases, as shown in Fig. [7] and Fig. [8] This 
phenomenon, again, may be understood based on the quasi¬ 
particle picture in Fig. [9] For two symmetric adjacent inter¬ 
vals, all the four kinds of entangled pairs contribute to £ for 
t < l. For two symmetric disjoint intervals, however, only the 
fast-fast pairs contribute to £ during d < t < d + l, and the 
other three kinds of entangled pairs do not make any contribu¬ 
tion at all. This explains why the numerical results are much 
smaller than the CFT results for the disjoint cases. 

The above quasi-particle picture on a lattice model may 
lead to many interesting results. For example, in the symmet¬ 
ric disjoint case, careful readers may wonder if the slow-slow 
pairs will make contributions to £ at a later time t > d + l. 
The answer is yes. Actually, as shown in Fig. [7] for the case 
d = 140 and l = 10, one can find that £ shows a “tail” or 
“revival” starting at around t = 180, which results from the 
contribution of slow-slow pairs. 

To demonstrate the above physical picture, we study, nu¬ 
merically, the distribution of entangled pairs as a function of 
(vl, vr) based on the harmonic chain. The strategy is as fol¬ 
lows. By choosing disjoint intervals A\ and A■> as shown in 
Fig J4] we fix the length of two intervals li = l 2 = l = 5 and 
time t = 100. Then we measure £ by changing the distances 
d\ and da separately. In this way, we plot £ as a function of 
{vl/v c ,vr/v c ) = (— di/t, da/t). As shown in Fig. 10 


it is 


found that there is apeak located at (—Vl/v c , vr/v c ) ~ (1,1) 
as expected from the CFT approach which is mainly con¬ 
tributed by the fast-fast pairs. On the other hand, there is an¬ 


other region ( — vl, vr < v c ) with finite £, which is mainly 
contributed by the other three kinds of entangled pairs, i.e., 
fast-slow pairs, slow-fast pairs and slow-slow pairs, respec¬ 
tively. While suppressed as compared to the main peak at 
(1,1), these quasiparticle pairs still carry non-negligible con¬ 
tributions to the entanglement negativity, leading to the tail or 
revival of the negativity £. 

Then a natural question is: how can we understand the 
“shape” of entangled pairs, i.e., the distribution of entangled 
pairs as a function of (y /■. v H ) ? In particular, for different 
methods of local quenches lITTl l22l 12814301 . the shape of en¬ 
tangled pairs may be very different. A quantitative study of 
this question as well as its application to the finite tempera¬ 
ture case will be reported elsewhere ED. 

Another question that the experts may ask is: For the time 
evolution of the entanglement negativity for two disjoint inter¬ 
vals after a global quench in Ref. ll26l . why cannot we observe 
an apparent mismatch between the CFT results and the nu¬ 
merical results? The reason is that Ref. ll26ll mainly focuses 
on the case l > d, which is close to the case of adjacent in¬ 
tervals. The shape of entangled pairs plays an important role 
only when d l, in which one can separate the contributions 
from the fast-fast pairs and the other three kinds of entangled 
pairs, as studied in our current work. To demonstrate this, we 
checked the entanglement negativity evolution of two sym¬ 
metric and disjoint intervals after a global quench in a har¬ 
monic chain. It is found that in the case d^> l, the numerical 
results are much smaller than the CFT results in magnitude, 
which agrees with our physical picture here. 
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IV. DISCUSSIONS AND CONCLUSIONS 


A. Relation between entanglement negativity and mutual 
information after a local quantum quench 


In the study of time evolution of the entanglement nega¬ 
tivity £{t) after a global quench l26l . it is found that £(t) 
shows the same feature as that of the Renyi mutual informa¬ 
tion / f r 'i (t). In our current study on the local quench problem, 
similar features are observed for the disjoint-interval cases. In 
this part, we hope to understand this observation, from a more 
general point of view, for both cases of adjacent and disjoint 
intervals after a local quantum quench. As a straightforward 
generalization of Ref. |26| . it is found that the relation be¬ 
tween the entanglement negativity and mutual information is 
independent of whether the quench of our interest is global 
or local. In other words, for both global and local quenches, 
the time evolution of the entanglement negativity £(t) has 
the same form as the Renyi mutual information /"(f) up to 
a global prefactor, as explicitly discussed in the following. 

The mutual information between two adjacent/disjoint in¬ 
tervals Ai and A 2 are defined as 


r(«) _ q(n) 1 o(") _ c(”) 

1 Ai,A2-—°Ai tD A 2 °A 1 uA 2 - 


(71) 


Based on the definition of Renyi entropy S^' 1 in Eq. (^j, we 
have 


constant. By neglecting various non-universal constants, then 
the Renyi mutual information in Eq. (73} can be expressed as 


r(«) 

1 A 1 ,A 3 


1 ln f ?7l,3772,4 \ 

71-1 \ 771 , 4772,3 / 


(76) 


By using the expression of A n in Eq. ( [T8| , one can immedi¬ 
ately get 


Al,A2 12 n v 771,3772,4/ 


(77) 


which has the same form as £ (f) in Eq. (42 1 apart from a con¬ 
stant prefactor |. Then the Renyi mutual information I^ua 2 
for two adjacent intervals can be obtained by taking the limit 
23 — f Z 2 - Since 


An) 

1 a 1 ua 2 


0 t ) 


2 Tl -|- 1 

3 n 


£AiUA 2 it) 


(78) 


for two disjoint intervals, one can get 


lim 

Z3->Z2 


tM 

A 1 U.A 2 


it) 


2n + l 
3 n 


lim £ AlU A 2 (t)- 

Z3~>Z 2 


(79) 


Although the above discussion is for the local quench, one 
can find that the same conclusion holds for the global quench 
ll26il . because our derivation is not sensitive to the concrete 
form of conformal mapping. 


/(") - = _1_ j n ( Tr/) TUA 2 \ p2) 

AlA ' n- 1 \ Trp n A Trp\ 2 ) ' ~ B. Conclusions 


Let us consider the disjoint-interval case first, i.e., A\ = 
[ 21 , 22 ] and A 2 = [ 23 , 2 : 4 ]. Expressed in terms of the corre¬ 
lation function of twist operators, we can get 


An) 

'Ai,A 2 


n — 1 


In 


(T n (21)7^(22)7^(23)7^(24)) 

(jn{zl)T n ( 22 )) (7^( 23 )^( 24 )) 


(73) 


By using the conformal mapping in Eq. (11 1 , we have 


(UrnMTM) 

i 

=n 

i 

where (• • • )rhp are correlators on the right half plane. The 
27V-point function of twist fields on the RHP reads 


dw 

dz 


Hi) 


<11 Tri( u 72i-l)7(i('U'2i))RHP 


(nr n Ki)r,K)W 

i 

C N n ( IIj<fc772fc,2j 772fc-l,2j-l 
n i\wi-wi\ A ™ l rij,fe772i-i,2fe 


A„ 

•^"({t lj,k})i 
(75) 


where the cross ratios rjij are defined in Eq. ( [26} . Here we are 
only interested in the limit rjij —>■ 0 or 1, i.eTT^ilrij^k}) is a 


In this paper, we studied the time evolution of the entangle¬ 
ment negativity that results from a local quench in conformal 
field theories, where the local quench is introduced by con¬ 
necting two decoupled CFTs. Once the two CFTs are joint 
at the endpoints, the interaction is simultaneously introduced, 
and then local quasiparticle excitations are generated at the 
jointing point. These quasiparticles carry information about 
entanglement and propagate freely along the system. The en¬ 
tanglement negativity of two intervals are built through these 
propagating quasiparticles. 

Then by employing CFT approach, we calculated the entan¬ 
glement negativity evolution for both adjacent intervals and 
disjoint intervals respectively. For the case of two adjacent 
intervals, the entanglement negativity grows as InT: initially, 
and then develops a plateau-like feature. When the quasiparti¬ 
cles propagate out of the intervals, the negativity drops to the 
ground-state value. For the case of two disjoint intervals, there 
is no entanglement negativity at the very beginning of local 
quench until the quasiparticles reach the two intervals simulta¬ 
neously. In the limit d l, a long-range entanglement which 
is independent of the distance d between two intervals is built 
through the quasiparticles. Then again, similar with the case 
of adjacent intervals, once the quasiparticles propagate out of 
the two intervals, there is a sharp drop of the entanglement 
negativity. 

Because our results are obtained from the CFT approach, 
the conclusion only applies to critical systems with a linear 
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dispersion relation in (1+1) dimensions. In other words, in 
the CFT approach, all the quasiparticles propagate at the same 
speed (“speed of light”). For a general lattice model such as 
the harmonic chain considered in this paper, however, the dis¬ 
persion relation is not linear for all momentum vectors. There 
are some quasiparticles propagating slowly compared with 
the speed of light. This is why the numerical results do not 
agree with the CFT results perfectly, although the main fea¬ 
tures agree with each other. 

Last but not least, we mention some interesting future prob¬ 
lems to be studied as follows. 

• Finite temperature effects. Most recently, finite tempera¬ 
ture effects on the entanglement negativity in conformal field 
theories were studied BT1 . In particular, finite temperature ef¬ 
fects on a local quench of the entanglement negativity were 
studied numerically based on a critical harmonic chain l27l . 
It is interesting to generalize our method to the finite temper¬ 
ature case, and have an analytical picture of the finite temper¬ 
ature effects on a local quench of the negativity. 

• Quantum entanglement of local operators, and entangle¬ 
ment density. In our work, a local quantum quench is real¬ 
ized by joining two separate CFTs at the endpoints. Another 
method to realize a local quantum quench is through acting 
with a local operator on an infinitely extended system [12814301 . 
By studying the increase of the Renyi entropy at a later time, 
quantum entanglement of local operators can be detected. In 
addition, by following the change of the entanglement entropy 
in a certain interval, one can study the entanglement density, 
which measures the density of entangled pairs between given 
two points E8l|42). Here, in our work, the local quench of 
the negativity provides a natural platform for studying the in¬ 
crease of entanglement, and therefore may be used to extract 
the entanglement of local operators, as well as the entangle¬ 
ment density. 

• Check lattice models. In our work, to demonstrate the 
CFT results, we study the critical harmonic chain numeri¬ 
cally. It will be interesting to check the CFT predictions in 
more complicated lattice models, such as itinerant fermions 
lfl8l and spin-chain systems, which are described by the Lut- 
tinger liquid theory. 

• Local quench of the entanglement negativity and entan¬ 
glement renormalization. Entanglement renormalization pro¬ 
vides a helpful framework to study the connection between 
quantum entanglement and its holographic geometry | [43l(46l . 
Global quenches and local quenches in entanglement renor¬ 
malization are discussed in Refs. B28l[45ll. respectively. In par¬ 
ticular, the effect of local quench on the entanglement entropy 
evolution based on entanglement renormalization is briefly 
discussed in Ref. (28]. It is of great interest to study the local 
quench of the entanglement negativity within the framework 
of entanglement renormalization. 

• Dynamical charged entanglement negativity Most re¬ 
cently, the concept of the charged Renyi entropies was pro¬ 
posed and studied in several works fiTi [48] . In particular, 
in Ref. (47] . the dynamical evolution of the charged Renyi 
entropies was studied. It is interesting to consider the entan¬ 
glement negativity in the presence of angular momentum and 
f/(l) charge, and study its dynamical properties after global 


quenches or local quenches. 

• Holographic study of time evolution of entanglement neg¬ 
ativity after a local quantum quench It is worth mentioning 
that the mutual information after a local quench has been stud¬ 
ied with holographic approach recently (28lf36l . On the other 
hand, holographic approach to the entanglement negativity 
was only discussed in the equilibrium case, although more 
needs to be understood (49l . It will be of great interest to study 
the non-equilibrium properties of the entanglement negativity 
with holographic methods. In addition, in a very recent work 
(50l on global quench, it is found that the free quasiparticle 
picture does not hold in a (1+1) dimensional CFT with c > 1 
(assuming no extended symmetry algebra). Relevant features 
were also observed in Ref. (361 . where it is found that the holo¬ 
graphic mutual information after a local quench shows qual¬ 
itatively different features from the CFT results. We believe 
that the same conclusion should hold in the local quench prob¬ 
lem of the entanglement negativity, and we leave this for our 
future work. 

Note added: After we finished this work, we noticed that 
related results on local quenches of the entanglement negativ¬ 
ity of two adjacent intervals appear in a recent workjST]. In 
Ref. EQ, the time evolution of the entanglement negativity 
after a local quench was studied in finite temperature, and the 
authors mainly focus on the adjacent-interval case. We hope 
our current work provides a platform for the study of disjoint- 
interval case after a local quench in finite temperature. 
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VI. APPENDIX 


In the appendix, we give explicit forms of Wij, which are 
used to calculate the entanglement negativity in the main text. 
The conformal mapping we used is 

where 


z = l + It. 


Then one has 


with 


l + IT 1 I-T. - Tl - —~ 

w = -1 —i/e 2 + (Z + it) 2 

1 


:=- (l + ir + pe l9 ) 


/9= [( e 2 + ;2_ T 2 ) 2 +4Z 2 r 2]l/4 ) 

1 2 It 

0 =-arctan—----. 

2 e 2 + l 2 - t 2 


(81) 


( 82 ) 
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The real time evolution can be obtained by replacing r with it 
in the last step. By expanding to the second order in e, p cos 9 
and psinZZ are expressed as follows. For l > 0, one has 


p cos 8 — »max[Z,t] 1 + 


2 [(max[Z,Z]) 2 — (min[Z,Z]) 2 


psin 8 —H min[Z, t] I 1 + 


2 [(min[Z,Z]) 2 — (max[Z,Z]) 2 


Based on the expressions of Wy, we can obtain the cross¬ 
ratio rjij as follows. 


912 - < t 2 e 2 


4 Z 2 (Z 2 -Z 2 )’ 


For l < 0, one has 


pcos 9 — > max[|Z|, t] I 1 + 


2 [(max[|Z| , t ]) 2 - (min[|Z|, Z]) 2 ] 


psm .0 — imin[|/|, t] I 1 + 


2 [(min[|Z|, Z]) 2 - (max[|Z|,Z]) 2 ] 


A. Appendix I: Wij for symmetric adjacent intervals 


9i3 = < re 2 
{ (Z 2 - Z 2 ) : 


-. t > l 


923 ~ < pe 2 


4 t 2 (t 2 — l 2 ) 


. t > l 


It is found that for the cases t <C l, t = l + 0 and t > l, one 
always has rpj = 1 or 0. 


B. Appendix II: w t:l for asymmetric adjacent intervals 


^(t + l) ’ 



-y/t(t — l), t > l 


— v l 2 — Z 2 , Z < 
e 


2Z V Z 2 — Z 2 ’ l> 


For asymmetric adjacent intervals, W \2 and w ^ have the 
same form as those in the symmetric case, and we only need 
to change l with l\. Similarly, for W 23 and w 2 3, we just change 
l with I 2 . Therefore, what we need to recalculate here are w \ 3 
and )/;,as follows. In the case of l\ < I 2 , one has 

-\Z?2 — Z 2 , t < 1 1 

2 

Wl3 = 7 V(h + h)(h ~ t), h<t<l 2 ( Q2 '1 


' (h+k ) 2 

(Z — l 2 )(Z + Zi) ’ 


= S 2 


W 23 = 


— v l 2 — Z 2 , Z < 
e 


— v Z 2 — Z 2 , Z > 
e 


-yl{l ~ t), t < l 


Z(Z — Z) ’ 


I — \J 1(1 + z), Z < l 

V 23 — j 2 

^ — \/z(z + z) . Z > l 


'll - z 2 , 


/ (Z + Z 2 ) {i 2 — Zi), l\ <t <l 2 
/ (z + Z 2 )(Z — Zi). Z > l 2 


In the case of Zi > l 2 , one has 

I {h + I2 ) (Z2 + Z) 
U, 13 = V (Z-Z 2 )(Zi+Z) ■ 

I (Zi + 1 2 ) 2 
. V (Z + Zi)(Z-Z 2 )’ 


12 < t < /l 
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'll ~ t' 2 , 


(h ~ h)(t + h) 
(h ~ t)(t — h) 


, l 2 <t <1 1 (95) 


~\/(t + h)(t — li). t > l\ 


-\/d 2 — t 2 , t < d 

e 


2d V )2Trf2> d<t<d + l 


Without loss of generality, here we choose l\ < h, as what 
we did in the main text. Then one has 


l\ — t 
h + t 1 
l 2 e 2 1 

4 't 2 {t 2 -i 2 y 
l 2 e 2 1 

4 't 2 {t 2 -l 2 Y 


l\ < t < l 2 


t 2 — d 2 


t d +1 


— \J{d + l + t)(d + l — t), 
2 

— \Jl(yd + l + i), 

2 

— \J{d + l + t)(t — d), 


d < t < d + l 


t d +1 


I 2 ~ t 
h + t ’ 
h ~ t 
h+t ’ 

1 

4 ' t 2 (t 2 -i 2 y 


l\ < t <L I 2 


— V (d + Z + t)(d -b Z — £), 
2 

~\/ (d + Z + t)(d + Z — t), 
2 

— ^/ (t + d + Z)(t — d — Z), 


d < t < d + Z 


t 'y> d + l 


ih-t){h+h) , . 1 

m 3 = < (h + t)(h-iiY 1 2 (98) 

e 2 (h+h) 2 .. 

(t 2 -i 2 ){t 2 -i 2 Y 

It is found that for the cases t -C Zi, t = h + 0“ and t > Z 2 > 
one always has rjij = 1 or 0. 


C. Appendix III: Wij for symmetric disjoint intervals 


■U’13 = 


2 ^/(d + t){d-t), 


l( t + d ){^r d -TTJTi)^ d<t< ‘ 
l {2d + l) {l^d-TVdTi)’ t> ' 


— y (d + Z + t) (d + Z — £), 
2 

— \/(d + Z + t)(d + Z — t), 


2 (d + Z) 4 


(* + d + Z)(t-d-Z)’ 


d < t < d + Z 


£ !> d -b Z 


-yd? — t 2 , 
€ 

-Vt 2 -d 2 , 

e 

-Vt 2 -d 2 , 

e 


d < t < d + l 


t ]> d ~b Z 


y/(d + t)(d-t), 

1 l(t + d) 

(t — d)(d + l — t) 


, d <t < d + l (100) 


- V (f + d)(t — d — l), t>d + l 


\i + d + t){i + d-t), 
/ (i + d-t){i + 2dy 


(2d + Z) * 


(t + d)(t — d—l) 


d <t < d + l 


t d + l 
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D. Appendix IV: Wij for asymmetric disjoint intervals 


<d\ - f 2 , 


-yd \ -1 2 , 


(di + d 2 + 0(^2 + t) 
(t — d 2 )(di + l + f) 

1 (dl+d 2 +Q2~ 

(£ + di + /)(£ — d 2 ) 

^ (dl + d 2 + £) 2 
(£ + di + Z)(t — d 2 ) 


t < di 
di < t < d 2 


d 2 < C £ <C di + l 


t > d 2 H - l 


\\Jd\- t 2 , 

- \/(di + d 2 )(d 2 — t), 

I (di + d 2 ) 2 ~ 

V (i - d 2 )(t + di) ’ 

I (di + d 2 ) 2 

V (i - d 2 )(t + di) ’ 

/ (di + d 2 ) 2 

V (* - d 2 )(t + di) ’ 


t < di 


di < t < d 2 


d 2 <C f < di + l 


t '> d 2 “t" l 


\\jd\ - 1 2 , 

^\/di -i 2 , 

I (di — d 2 + ()(£ “1“ d 2 ) 
Y (di + l — t)(i - d 2 ) 

- \/(i + d 2 )(t — di — I), 
-\/{t + d 2 )(t — di — Z), 


t < di 
di < t < d 2 




t > d 2 “t - l 


-y d 2~ t2 ’ 


t < di 


^\/{t + d 2 )(d 2 - di), d 1 <t<d 2 

- \/(t + d 2 )(t - di), d 2 < t < di + l 

- £ \/(t + d 2 )(t-di), di + l < t < d 2 + l 

2 

- \/(t + d 2 )(t - di), t > d 2 + l 

( 112 ) 


/ (d 2 + /) 2 -f 2 , 
/ (d 2 + /) 2 -f 2 , 
/ (d 2 + /) 2 -f 2 , 


- V (di + d 2 “h ‘2l)(d 2 + l — £), di 

I (di + d 2 + 2/) 2 
V (t — cZ 2 — 0(^ + + 0 ’ 


£ < di 
di < £ < d 2 


di + / 


£ do + / 


- V (d2 + Z) 2 — £ 2 , £ < di 

- a/( di + c^2 H - i){d 2 + Z — £), di < £ < c?2 

2 

- \/(di +6^2 + l)(d 2 + Z — £), d2 < £ < di + / 
2 

-\/(di + d2 + l){d2 + l — £), di + / < t < d2 + / 


(di + d2 + /) 2 
(t- d 2 - l)(t + di) ’ 


t d 2 + / 


— v(d 2 + 0 2 — 

^\/(d 2 + l) 2 - t 2 , 

- £ V(d 2 + l) 2 ~ti, 

2 

— \/(t + d 2 + ?) (d 2 — di), di 

2 

— \/(f + d 2 + ?)(i — di — /), 


t < di 


d\ < t < d 2 


d 2 <C t <C di + / ^24 — 


t > d 2 + / 
(HO) 


— y(d 2 + l) 2 — t 2 , t < di 

2 

— ^/(t + d 2 + l)(d 2 + Z — di), di < t < d 2 

— a/( t H - d 2 + l)(d 2 + Z — di), d 2 < t < di + l 
2 

— a/(£ H - d 2 + l)(d 2 + Z — di), di + / < t < d2 + / 
2 

— a/(^ H - d 2 + /)(t — di). t > d 2 + l 

( 114 ) 
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